Comment on "Critique of g-entropy for thermal statistics" by M. Nauenberg 



cn 

O 

o 

Oh: 

<^ 

o 



I 

o 
o 



> 

O 
m 
o 



I 

o 
o 



X 



C. Tsallis * 

Centra Brasileiro de Pesquisas Fisicas, Rua Xavier Sigaud 150, 22290-180 Rio de Janeiro - RJ, Brazil 

(February 2, 2008) 

Pacs Numbers: 05.70.-a, 05.20.-y, 05.90.+m 



It was recently published [1] a quite long list of ob- 
jections about the physical validity for thermal statistics 
of the theory sometimes referred to in the literature as 
nonextensive statistical mechanics. This generalization 
of Boltzmann-Gibbs (BG) statistical mechanics is based 
on the following expression for the entropy: 



S = k ~^^fP- (q eTZ;Si = Sbg = -fc In p.) . 

(1) 

The author of [1] already presented orally the essence 
of his arguments in 1993 during a scientific meeting [2]. 
I am replying now simultaneously to the just cited pa- 
per, as well as to the 1993 objections (essentially, the 
violation of "fundamental thermodynamic concepts" , as 
stated in the Abstract of [1]). The list of objections and 
critical statements being extremely long, it is perhaps 
not really necessary at the present stage to reply to all 
the points. For time and space economy, I will therefore 
address here only a few selected points, hopefully the 
most relevant ones, physically and/or logically speaking. 

About the nonextensivity of the entropy Sqt 

The entropy Sq is nonextensive for independent sys- 
tems (see Eq. (6) of [1]), which by no means implies 
that it cannot he extensive in the presence of correlations 
at all scales. Nowhere in [1] is there clear evidence of 
taking this fact into account in what concerns the valid- 
ity of the g-thermostatistics. It is nevertheless of crucial 
importance, as we illustrate now for the simple case of 
equiprobability (i.e., pi = l/W^ Vi). In such simple situ- 
ation, Eq. (1) becomes 



w 



Sg ^ klllqW 



{lnqX = (x^ ' — 1)/(1 — (7); Ini a; = Inx) . 

(2) 



If a system constituted by N elements is such that it can 
be divided into two or more essentially independent sub- 
systems (e.g., independent coins or dices, or spins inter- 
acting through short-range coupling), we generically have 
W ^ fi^ {fi> 1). Consequently, Sq/k - ln,/i^. There 



is an unique value of q, namely q — 1, for which we obtain 
the usual result Sq oc N. But if the system is such that 
we have W NP {p > 0), then Sq/k ~ In, TV. Once 
again, there is a unique value of g, namely q = 1 — 1/p 
for which, 5*^ oc N. This fact is well known to many 
scientists working on nonextensive statistical mechanics, 
and has been published in the special volume dedicated 
to the subject indicated in Refs. [5, 6, 14] of paper [1]. 
The same property holds in fact for Sc,(^q) , a-{q) being any 
smooth function of q such that cr(l) = 1 (e.g., a = 1/q, 
OT a = 2 — q). For the correlated case, we have S„(^q) oc N 
only for q satisfying {1 — a{q))p = 1. The relevance of 
this property (S oc N) for thermodynamics needs, we 
believe, no further comments. 

About the concept of "weak coupling" in [I]: 

Much of the criticism in [1] involves the concept of 
"weak coupling" . To make this point clear through an 
illustration, let us think of the ground state of a Hamil- 
tonian many-body classical system whose elements are 
localized on a d-dimensional lattice and have two-body 
interactions among them. Let us further assume that the 
(attractive) coupling constant is given by Cy = — c/r^- 
(c > 0, a > 0, and = 1,...). The potential en- 
ergy U{N) per particle generically satisfies U{N)/N oc 



-c dr r"^' 



Therefore, for a/d > 1 [short-range interactions in the 
present context), we have that limf^^ooU{N)/N is finite, 
and BG statistical mechanics certainly provides the ap- 
propriate answer for the stationary state (thermal equi- 
librium) of the system. In this case, all the usual pre- 
scriptions of thermodynamics are satisfied, as well known 
[3]. If the interactions are, however, long-range (i.e., 
< a/d < 1), then limN^ooU{N)/N diverges, and the 
case needs further discussion. It might well happen that, 
dynamically speaking, the N ^ 00 and the < — > 00 limits 
do not commute. If so, only the limjv^oo linit^oo or- 
dering corresponds to the BG stationary state, whereas 
the opposite ordering, limt^oo liniAr^ooj might be a com- 
plex one, different from the BG state, and in some oc- 
casions possibly related to the one obtained within the 
g- formalism. It is clear then that, if we have long-range 
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interactions and >> 1 (say of the order of the Avo- 
gadro number), it might very weU happen that the BG 
equiUbrium is physically inaccessible, and the only physi- 
cally relevant stationary or quasi-stationary (metastable) 
state is a non-Gibbsian one. Such situation is indeed 
found in [4], as discussed below. 

We can now address the manner used in [1] to refer 
to "weak coupling" . It applies essentially in the simple 
manner stated in [1] only for a/d > 1, being conceptually 
much more subtle for < a/d < 1. For example, if 
< a/d < 1, U{N)/N diverges as N^-"/^ {N oo) for 
any nonvanishing value of c, even for c corresponding to 
... lO^^^eVl Consistently, the generic use, without fur- 
ther considerations (such as the {N,t) (oo,cx)) limits, 
and the range of a/d), of relations such as Eqs. (5) and 
(7) of [1] seems irreducibly unjustified; as they stand, 
they trivially yield to no other possibility than q = 1. In 
fact, this point has already been transparently addressed 
by Fermi in 1936 [5]. 

About the determination of the value of q for a given 
system: 

The entropic parameter q is referred in [1] as an "unde- 
termined parameter" . Moreover, the author claims hav- 
ing proved that "g must be a universal constant, just 
like the Boltzmann constant A;...". I have difficulty in 
unambiguously finding in the paper whether this kind of 
statement would only apply to Hamiltonian systems, or 
perhaps also to dissipative ones; to systems whose phase 
space is high-dimensional, or perhaps also to the low- 
dimensional ones. By "undetermined" , it remains not 
totally clear whether the expression is used in the sense 
that q is "undeterminable" , or in the sense of "not yet 
determined" . However, if we put all this together, one 
might suspect that what is claimed in [1] is that it can 
be determined from first principles, and that the author 
has determined it to necessarily he q — I. 

To make this point transparent, we may illustrate the 
factual nonuniversality of q by addressing the logistic- 
like family of maps Xt+i = 1 — a|x(|^, whose usefulness 
in physics can hardly be contested (at least for z = 2). 
As conjectured since 1997 [6], numerically exhibited in 
many occasions (e.g., in [7-10]), and analytically proved 
recently [11,12] on renormalization group grounds, q does 
depend on and is therefore not universal, in neat con- 
trast with what is claimed in [1]. Its value for z = 2 (i.e., 
the standard logistic map) , as given by the sensitivity to 
the initial conditions, is q — 0.244487... at the edge of 
chaos (e.g., a = 1.401155...), whereas it is g = 1 for all 
values of a for which the Lyapunov exponent is positive 
(e.g., for a = 2). We have illustrated the nonuniversal- 
ity of q for nonlinear dynamical systems with its value 
at the edge of chaos of the logistic map. It is perhaps 
worthy to notice that, since it has been proved to be an- 
alytically related to the Feigenbaum universal constant 
ap [1/(1 — (?) = lnaF/ln2], and since this constant is al- 
ready known with not less than 1018 digits, we actually 



know this particular value of q with the same number of 
digits. Such a precision is self-explanatory with regard 
to the fact that q can be determined from first principles 
and that it can be different from unity (see also [13,14]). 

A second illustration of the nonuniversality of q can be 
found in the three-component Lotka-Volterra model in a 
d- dimensional hypercubic lattice [15]. This illustration 
is quite interesting because this model is a many-body 
problem. The corresponding growth of droplets has been 
shown to yield, through imposition of the finiteness of 
the entropy production per unit time, q = 1 — 1/d for 
d = 1,2 [16]. This law has also been checked for d = 3, 4 
[17]. It is of a rather simple nature, essentially related to 
the fact that the growth of the bulk regions of this specific 
model is characterized by the droplet linear size linearly 
increasing with time [17]. A similar law is obtained for a 
quite different model, namely a Boltzmann d-dimensional 
lattice model for the incompressible fluid Navier-Stokes 
equations. Indeed, an unique entropy, namely Sq, with 
an unique value of q, is mandated by the imposition of 
the most basic Galilean invariance of the equations. For 
the single-speed single-mass model it is = 1 — 2/d [18]. 
For more sophisticated models q is determined by a tran- 
scendental equation [19]. These examples show how the 
entropy Sq enables to put on equal grounds situations 
that are physically quite disparate. 

Given the preceding illustrations of dissipative sys- 
tems, and many others existing in the literature, it could 
hardly be a big surprise if, also for many-body Hamil- 
tonian systems, q turned out to be a nonuniversal in- 
dex essentially characterizing what we may consider as 
nonextensivity universality classes (in total analogy with 
the universality classes that emerge in the theory of crit- 
ical phenomena). More precisely, one expects q = 1 
for short-range interactions {a/d > 1 in the example we 
used earher), and q depending on (d, a) (perhaps only on 
a/d) for long-range interactions (i.e., < a/d < 1), in 
the physically most important ordering limf^oo liniAr— >oo- 
Although expected, the uncontestable evidence has not 
yet been provided. It is not hard for the reader to imag- 
ine the analytic and computational difficulties that are 
involved. However, suggestive results are accumulating 
which point towards the applicability of nonextensive sta- 
tistical mechanics for such long-range Hamiltonian sys- 
tems. Although we shall later come back onto this prob- 
lem, let us already mention the following points. 

(i) The one-body marginal distribution of velocities 
during the well known longstanding quasi-stationary 
(metastable) state of the isolated classical inertial 
XY ferromagnetically coupled rotators localized on 
a d-dimensional lattice can be anomalous (i.e., non- 
Maxwellian). Indeed, it approaches, for a non-zero- 
measurc class of initial conditions of the a — (yd) 
model and not too high velocities, a g-exponential dis- 
tribution (we remind that = [1 -|- (1 — q)x]^/^^~'^\ 
hence = e^) with (j > 1 [4]. If the energy dis- 
tribution followed BG statistics, the one-body marginal 
distribution of velocities ought to be quasi-Maxwellian 
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(strictly Maxwellian in the N ^ oo limit since then the 
microcanonical-ensemble necessary cutoff in velocities di- 
verges), but it is not . As specifically discussed in [4], 
the numerical results are incompatible with BG statis- 
tics. However, they do not yet prove that the one-body 
distribution of velocities precisely is, for the canonical en- 
semble, the one predicted by noncxtensivc statistics. In- 
deed, considering the appropriate limit (iV, M, N/M) 
(cxD, oo, oo) {N being the number of rotators of the iso- 
lated system, and M being that of a relatively small sub- 
system of it) is crucial. Work along this line is in progress. 

(ii) In the same model, at high total energy, the 
largest Lyapunov exponent vanishes like where k 
depends on a/d [20,21]. Also during the lonstanding 
state, the largest Lyapunov exponent vanishes, this time 
like [22]. It is clear that, with a vanishing Lya- 
punov spectrum, the system will be seriously prevented 
from satisfying Boltzmann's "molecular chaos hypothe- 
sis", hence the "equal probability" occupation of phase 
space. 

(iii) In the longstanding regime of the a = (Vd) 
model, there is aging [23] , something which is totally in- 
compatible with the usual notion of thermal equilibrium. 
The correlation functions depend on the "waiting time" , 
and are in all cases given by g-exponential functions. 
Even at high total energy, where the one-body distri- 
bution of velocities is Maxwellian, and where there is no 
aging, the time correlation functions are still given by q- 
exponentials with q > 1, instead of exponentials, which 
is the standard expectation in BG statistics. 

(iv) The temperature (oc mean kinetic energy per par- 
ticle) relaxes, after the quasistationary state observed in 
the one-dimensional < a < 1 model, towards the BG 
temperature through a g-exponential function with q> 1 
[24]. 

(v) In Lennard- Jones clusters of up to A'^ = 14 atoms, 
the distribution of the number of links per site has been 
numerically computed [25], where two local minima of 
the many-body potential energy are "linked" if and only 
if they are separated by no more than one saddle-point. 
This distribution is a g-exponential with g ~ 2, as can be 
checked through direct fitting. The possible connection 
with our present discussion comes from the fact that the 
average diameter of the cluster is (in units of atomic size) 
of the order of 14^/'^ ~ 2.4. Consequently, although the 
Lennard-Jones interaction is not a long-range one ther- 
modynamically speaking (indeed, a/d = 6/3 = 2 > 1), 
it can effectively be considered as such for small clus- 
ters, since all the atoms substantially interact with all 
the others. 

(vi) The distribution of the number of links per node 
for the Albert-Barabasi growth model [26] yielding scale- 
free networks is analytically established to be, in the sta- 
tionary state, a g-exponential with g = [2m(2 — r) + 1 — 
p — r]/[m{3 — 2r) + l—p — r] > 1, where {m,p, r) are mi- 
croscopic parameters of the model. If we associate to this 
network an a = interaction per link, the just mentioned 
distribution also represents the distribution of energies 



per node. Although this is not the same distribution as 
that of the energy of microscopic states associated with 
a Hamiltonian, it is neither very far from it. 

(vii) Although not being many-body problems, let us 
mention at this point some results that have been ob- 
tained with the d = 2 standard map and with a d = 4 
set of two coupled standard maps. Both systems are 
conservative and simplectic, having therefore the dynam- 
ical setup of a standard Hamiltonian. The d = A sys- 
tem has Arnold diffusion as soon as the nonlinear cou- 
pling constant a is different from zero; this guarantees a 
chaotic sea which is singly connected in phase space (we 
may say that Oc = 0). The structure is more complex 
for the d = 2 case because no such diffusion is present; 
consistently, unless a is sufficiently large, disconnected 
chaotic "lakes" are present in the phase space; below 
a.c = 0.97..., closed KAM regions emerge in the prob- 
lem. The remark that we wish to do here is that, in 
strong analogy with the many-body long-range Hamilto- 
nian cases we have been discussing, both the d = 2 and 
the d = 4 maps present a longstanding quasistationary 
states before crossing over to the stationary ones. The 
crossover time tcrossover diverges when a approaches Oc 
from above. This is very similar to what happens with 
the above {d, a) Hamiltonian, for which strong numeri- 
cal evidence exists [4,22,28] suggesting that tcrossover di- 
verges as (N^-"/'^ - 1)/(1 - a/d) when N —>■ oo. 

Although none of the (seven) factual arguments that 
we have just presented constitutes a proof, the set of 
them does provide, in our undestanding, a quite strong 
suggestion that the longstanding quasistationary states 
existing in long-range many-body Hamiltonians might 
be intimately connected to the nonextensive statistics, 
with g depending on basic model parameters such as d 
and a. The entropic index g would then characterize 
universality classes of nonextensivity, the most famous 
of them being naturally the g = 1, extensive, universality 
class. Such viewpoint is also consistent with the discus- 
sion about non-Gibbsian statistics presented in [29]. Last 
but by no means least, it is consistent with Einstein's 
1910 criticism [30] of the Boltzmann principle S = klnW 
(lengthily commented in Ref. [6] of [1]). 

About thermal contact between systems with different val- 
ues of q and the 0*^ principle of thermodynamics: 

We focus now on a strong and crucial statement in 
[1], namely "... a Boltzmann- Gibbs thermometer would 
not be able to measure the temperature of a g-entropic 
system, and the laws of thermodynamics would therefore 
fail to have general validity." [1]. We shall present here 
the results [39] of molecular- dynamical simulations (us- 
ing only F = ma as microscopic dynamics) which will 
precisely exhibit what is claimed in [1] to be impossible. 
We shall illustrate this with the isolated a = model of 
planar rotators, and proceed through two steps. 

We first show (Fig. 1) how the "temperature" (defined 
as twice the instantaneous kinetic energy per particle) of 
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a relatively small part of a large system relaxes onto the 
"temperature" of the large system while this is in the qua- 
sistationary regime (where the system has been definitely 
shown to be non-Boltzmannian, and where it might well 
be described by the g-statistics) . We verify that the rest 
of the system acts for a generic small part of itself as a 
"thermostat" , in total analogy with what happens in BG 
thermal equilibrium. This is quite remarkable if we think 
that tli(' system is in a state so different from thermal 
equilibrium that it even has aging! 

We then show (Fig. 2) how a BG thermometer (its 
internal degrees of freedom are those of first-neighbor- 
coupled inertial rotators, hence definitively a q ~ I sys- 
tem) does measure the "temperature" of the infinitely- 
range-coupled inertial rotators during their quasistation- 
ary state, hence where the statistics is definitely non- 
Boltzmannian. At the light of this evidence, it appears 
that the 0*'' principle of thermodynamics is even more 
general than the already important role that BG statisti- 
cal mechanics reserves for it. Naturally, the fluctuations 
that we observe in both figures are expected to disappear 
in the {N, M, N/M) —>■ (oo, oo, oo) limit. 

The facts that we have mentioned up to this point 
heavily disqualify the essence of the critique presented 
in [1]. I believe, nevertheless, that it is instructive to 
further analyze it. 

About the existing mathematical foundations of nonex- 
tensive statistical mechanics: 

It is essentially claimed in [1] that it can be proved, 
from the very foundations of statistical mechanics, that 
the only physically admissible one is that of BG. It is 
however intriguing how such a strong statement may be 
done without clearly pointing the mathematical errors 
that should then exist in the available proofs of the q- 
exponential distribution. Such proofs have been provided 
by Abe and Rajagopal [31-34]; they are multiple, mu- 
tually consistent; and generalize the well known proofs 
done, for BG statistics, by Darwin- Fowler (in 1922), 
Khinchin (in 1949) and Balian-Balazs (in 1987), respec- 
tively using the steepest-descent method [31], the laws of 
large numbers [32] , and the counting for the microcanoni- 
cal ensemble [33]. All these proofs are ignored in [1]. The 
critique therein developed outcomes severely diminished. 

Similarly, no mention at all is made in [1] of the 
^-generalizations of Shannon 1948 theorem, and of 
Khinchin 1953 theorem, which are universally considered 
as part of the foundations of BG statistical mechanics 
since they prove under what conditions Sbg is unique. 
These two ^-generalizations [35,36] analogously exhibit 
the necessary and sufficient conditions associated with 
the uniqueness of Sq. 

Finally, no mention at all is made of the fact that Sg 
(Vg > 0) shares with Sbg three remarkable mathematical 
properties that arc quite hard to satisfy, especially simul- 
taneously. These three properties are concavity (Ref. [1] 
of [1]), stability [37], and finiteness of entropy production 



per unit time (sec [38], among others). The difficulty of 
having such agreable mathematical features can be mea- 
sured by the fact that Renyi entropy (Eq. (19) of [1]), 
for instance, satisfies none of them for abitrary q > 0. 

It is perhaps for not paying due attention to all these 
theorems that the cyclic argument involving Eqs. (22- 
26) of [1] has been included in the critique. Indeed, 
that argument uses Eq. (22) to "prove" Eq. (26). Such 
a consistency can hardly be considered as surprising 
since the distribution in Eq. (22) is currently estab- 
lished precisely using the BG entropy, i.e., the form of 
Eq. (26). By the way, immediately after Eq. (26) we 
read "provided /(I) = /(O) = 0, which corresponds to 
the requirement that the entropy vanishes at T = 0". 
It is in fact only /(I) = which is related to the van- 
ishing entropy at T = 0. The property /(O) = has 
in general nothing to do with it; it is instead related 
to the expansibility of the entropy, i.e., the fact that 
S{pi,p2, ...,pw,0) = S{pi,p2, ■■■,pw)- 

About existing exact solutions of anomalous Fokker- 
Planck and Langevin equations: 

The standard = 1 Langevin equation (with a drift 
coefficient 7 and an additive noise), and the standard 
d = 1 Fokker-Planck equation admit as exact solutions 
the Gaussian distribution, and are usually considered 
as paradigmatic mesoscopic descriptions associated with 
BG statistical mechanics. They can be naturally gen- 
eralized by also including a multiplicative noise (with 
amplitude M) in the Langevin equation, and by con- 
sidering the so called "porous medium equation", i.e., a 
nonlinear Fokker-Planck equation where the Laplacian 
operator applies to the power of the distribution. The 
exact solutions of these two nontrivial (nonlinear) equa- 
tions are g-Gaussians, with q = (3il/ + 7)/(Af -I- 7) > 1 
for the former [40], and g = 2 — 1/ < 3 for the latter [41]. 
These suggestive mathematical facts are ignored in [1]. 

Miscellanous 

The precise formulation of nonextensivc statistical me- 
chanics has, since 1988, evolved along time in what con- 
cerns the way of imposing the auxiliary constraints un- 
der which Sq is optimized (see Refs. [1-3] of [1]). The 
paradigmatic case occurs for the canonical ensemble, 
where one must decide how to generalize the traditional 
energy constraint. The correct manner is nowadays ac- 
cepted to be that indicated in Ref. [3] of [1], i.e., Eq. (3) 
of [1], namely 



U„ 



(3) 



This particular writing of the energy constraint has var- 
ious interesting features. Let us mention here three of 
them (further convenient features can be found in Ref. 
[3] of [1]). 
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(i) It is precisely this form which emerges naturally 
within the steepest-descent proof [31] of the g-statistics. 
It is a trivial consequence of the fact that de^ /dx = (e^ . 

(ii) This particular form makes the theory to be, in 
what concerns the energy distribution, valid up to a sin- 
gle value of q, namely precisely that determined by the 
trivial constraint Yl^iPi ~ ^- illustrate this in 
the continuum limit, for a typical example where the 
density of states g{e) oc e'' for e — > oo (7 G 7?.). Since 
we wish p{e) to be normalizable, we must impose that 
/r„.ta„t*5(e)p(e). Since p{e) oc eV(i-9) for e ^ 00, it 
must be 7 + 1/(1 — q) > —1, hence 

g< (2 + 7)7(1 + 7) , (4) 

(g < 2 for the simple case of an asymptotically constant 
density of states, i.e., 7 = 0). The finiteness of con- 
straint (3) imposes /c'^„sto„t <^efi'(e) e b(^)]^ to be finite, 
which, interestingly enough, yields the same upper bound 
as before, namely Eq. (4). In other words, this makes the 
theory to have both constraints (norm and energy) math- 
ematically well defined (i.e., given by finite numbers) all 
the way up to a single upper bound for q. 

(iii) This structure (based on escort distributions [42]) 
for the energy constraint allows the construction of a 
quite general entropic form [43] which is extremized by 
the Beck-Cohen superstatistics [44], and which, quite 
remarkably, is stable [45] (like Sq, and in variance with 
Renyi entropy). 

Since 1988, many applications have been proposed for 
the nonextensive statistics. Some of these have been 
elaborated within the 1988 way of writing the energy 
constraint (Ref. [1] of [1]), others have been elaborated 
within the 1991 way of writing this constraint (Ref. [2] 
of [1]), and finally others with the 1998 way (Ref. [3] of 
[1]). It is unfortunate that the 1998 way was not found 
from the very beginning in 1988, but this is the way it 
did happen (for a variety of reasons that are essentially 
commented in Ref. [3] of [1]). Consistently, it seems fair 
to nowadays restrict possible criticism to applications in- 
deed using the 1998 version. Unfortunately, all types of 
applications are criticized in [1] independently from what 
particular manner have the author (s) adopted for the en- 
ergy constraint. An intriguing example of such procedure 
is the criticism of some 1995-1996 papers (Refs. [19-21] 
of [1]) on the possible ^-generalization of the black-body 
radiation law. They indeed satisfy, as they should and as 
claimed in [1], the Stefan-Boltzmann law {explicitly 
written in Eq. (18) of Ref. [19] of [1]). Nevertheless, 
they do not escape the criticism! It is argued in [1] that 
errors have been done in these three papers, and that, if 
these errors had not been done, the papers would have 
violated the T"^ law, and therefore they also deserve criti- 
cism. In addition to this somewhat courageous comment, 
no remarks are done about the fact that all three were 
published up to three years before the need for re- writing 
the energy constraint became clear. More significantly. 



this specific criticism is indeed intriguing since it can be 
trivially shown that the proportionality law remains 
the same for all energy statistical distributions (hence not 
only the EG one) as long as the microscopic energy scales 
linearly with the temperature (i.e., for photons, as long 
as the distribution depends on the light frequency i' and 
the appropriate temperature T, only through p/T). Only 
the proportionality coefiicient of the law depends on 
the specific statistics. 

The author of [1] claims to have delivered the epis- 
temological coup de grace to nonextensive statistical 
mechanics. Indeed, expressions like "unphysical" , "man- 
ifestly incorrect", "devoid of any physical meaning", 
"do not have any physical meaning" , "disregarding such 
basic considerations" , "nonsensical" , "failure of this for- 
malism" , "inconsistencies" , "inconsistent with the fun- 
damental principles of thermodynamics and statistical 
mechanics", "absolutely no physical justification has 
been given", and analogous ones, have been profusely 
used in [Ij. We have essentially argued here that what 
we are facing is rather the opposite, in the sense that it is 
precisely the basis of the critique in [1] which appears to 
be deeply inconsistent with very many, and by now well 
established, physical and mathematical facts. We have 
only addressed the main mispaths and inadvertences 
in [1]. There are several more, but the full consider- 
ation of them all would demand an appreciable effort 
which, at the present moment, does not seem worthy. 
Our overall conclusion is that, although several impor- 
tant and/or interesting points related to nonextensive 
statistical mechanics still need further clarification, this 
theory undoubtedly exhibits nowadays a sensible number 
of physically and mathematically consistent results. Of 
course, as it has always been, only time will establish 
its degree of scientific utility in theoretical physics and 
elsewhere. 
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FIG. 1. Time evolution of the temperature Tmicrocanonicai 
= 2K{N)/N {K{N) = total kinetic energy) of one isolated 
system started with waterbag initial conditions at (conve- 
niently scaled) energy per particle equal to 0.69 (A^ — 5000 
rotators; green line), and of the temperature T^anonicai = 
2K{M)/M {K(M) = subsystem total kinetic energy) of a 
part of it (M = 500 rotators; blue line) . The M rotators were 
chosen such that their temperature Tcanonicai was initially 
below (a) or above (b) that of the whole system. It is partic- 
ularly interesting the fact that, in case (b), the temperature 
of the subsystem of M rotators crosses the BG temperature 
Tbg = 0.476 without any particular detection of it. 



FIG. 2. Time evolution of the temperature Tthermoatat = 
2K{N)/N {K{N) = thermostat total kinetic energy) of 
one infinitely-range-coupled large system (thermostat) started 
with waterbag initial conditions (A'' = 100000 rotators; green 
line) and of the temperature Tthermometer = 2K(M)/M 
{K{M) = thermometer total kinetic energy) of one 
first-neighbor-coupled relatively small system (thermometer) 
started at Maxwellian equilibrium at a temperature below 
that of the thermostat (Af — 50 rotators; blue and red lines). 
The large system is in the quasistationary state (where it is ag- 
ing!); its (conveniently scaled) energy per particle equals 0.69. 
The thermometer-thermostat contact is assured by only one 
bond per thermometer rotator, and starts at time tcontact- 
The intra-thermostat and intra-thermometer coupling con- 
stants equal unity; the thermostat-thermometer coupling con- 
stant equals 0.001. The thermalization of the thermome- 
ter occurs at the thermostat temperature, and up to time 
t = 3 X 10^, exhibits no detection of the BG equilibrium 
temperature Tbg = 0.476 . The same phenomenon with the 
thermometer initial temperature being larger than that of the 
thermostat is not shown, because our numerical results sug- 
gest that the >> Af >> 1 limit has to be satisfied in an 
even more stringent manner due to the relatively large fluc- 
tuations of Tthermometer ■ FoT clarity, uot all the points of 
the curves have been represented, but they have been instead 
logarithmically decimated. 
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